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Abstract. We study various mathematical aspects of discrete models on graphs, specifi- 
cally the Dimer and the Ising models. We focus on proving gluing formulas for individual 
summands of the partition function. We also obtain partial results regarding conjectured 
limits realized by fermions in rational conformal field theories. 



1. Introduction 

The idea that the dimer and Ising statistical models on finite graphs have, as a limit, the 
free fermion conformal field theory in dimension 2 and 1, respectively, has become a well 
established theme in mathematical physics (cf. [29], [8]). Yet, great challenges remain on 
the road to developing this into a mathematically rigorous theory. A part of the difficulty 
is in reconciling the languages and concepts of graph theory and conformal field theory: the 
structure of conformal field theory predicts many features which are not readily visible on 
graphs. One must develop mechanisms how such features originate on the graph level and 
will fully manifest themselves in the limit. The purpose of this paper is, in some sense, 
to begin this investigation in earnest. We record some of the phenomena which need to 
modelled. We also obtain concrete results as first steps in the desired direction. 

While conformal field theory was long considered a part of the domain of physics, rigorous 
mathematical approaches now do exist. Perhaps most appealing is the "naive" approach 
due to Graeme Segal ^3]J, which asks for a Hilbert space H, and for each Riemann surface 
S with analytically parametrized boundary components, an element UiTi) (defined up to 
scalar multiple) of the Hilbert tensor product of copies of H and its dual indexed by the 
boundary components (depending on orientation), satisfying appropriate "gluing axioms". 
The paper [31] is just an outline, and details have been since filled in, in part by Kriz and 
his co-authors [2211181 [TT]. 



A competing approach, starting with the notion of vertex algebra, was started by Borcherds 
[3] and Frenkel-Lepowsky-Meurman |13] (see also Beilinson-Drinfeld [2j). This approach 
models only a part of the gluing structure postulated by Segal, but has the advantage that a 
vertex algebra is a completely algebraic concept, which does not require analysis. A vertex 
algebra, however, models only a part of conformal field theory, the "chiral 0-sector". An 
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extension of the approach which models a whole conformal field theory, was developed by 
Huang and Lepowsky [T9j, but it does involve analysis, although not Hilbert spaces. Efforts 
to model the entire conformal field theory as a completely algebraic concept are underway 

[171123]. 

The mathematically richest case of conformal field theory is a rational conformal field 
theory, when there is a finite set of labels A with an involution (?)* and the Hilbert space 
breaks up as 

(1) H = ^H^mx^. 

aga 

The Hilbert spaces Hx are called chiral sectors and are subject to additional axioms. Parallel 
to the formula ([T]), there is a formula on graded dimensions, the partition function formula 
(formula fHSj) below). A major advantage of the vertex algebra approach is that it predicts 
from a vertex algebra alone (i.e. algebraic data) when we get a rational conformal field 
theory, and gives an algebraic formula for its sectors (see also [2J). 

The free fermion is now completely described mathematically as a conformal field theory 
both in the Segal and vertex algebra formalisms (see [211 [221 US])- Yet, the theory is more 
subtle than one might think. We describe some of the issues involved in Section jHl below. 
For one thing, a rational conformal field theory in the precise sense outlined above actually 
describes the bosonisation of the free fermion. There is a variant of the above description, 
called a fermionic rational conformal field theory, which is more closely related to the topic 
of this paper. It is actually the fermionic rational conformal field theory analogue of the 
partition function formula (formula ( 15^ . and its analogues for higher genus) we are interested 
in. In the case of the free fermion, this formula (fTOj) is precisely a sum over spin structures 
of a Riemann surface of genus g, with a factor 1/2^, which matches the formula on graph 
dimer partition function found by Cimasoni and Reshetikhin [Ij. 

The present paper is mainly about graphs. How can we model the concepts mentioned 
discretely? Fortunately, great strides have already been made. Mercat and Kenyon (201 [2T| 
defined critically embedded graphs on Riemann surfaces, which give discrete models the 
conformal structure. A basic quantum field theory-like gluing property of the dimer model 
was observed by Cimasoni and Reshetikhyn [3]. 

A large part of the present paper (Sections [21 [3|) consists of examining the gluing property 
in more detail, both for the Ising and dimer models. Our main new result in this direction is 
proving a gluing formula not just for the entire partition function, but for its summands, i.e. 
discrete analogues of ( ITO|) . (15^ . More precisely, we consider arbitrary labellings, not only 
critical ones; this generality, the formula is actually related to quantum field theory (not 
conformal field theory). The conformal case, however, is the one we need to consider for 
taking scaling limits. For simplicity, we restrict to the case of a double torus being cut into 
two surfaces of genus 1, but it is straightforward to generalize the formula to more general 
cuts along separating curves. 

Our gluing formula is substantially more difficult to obtain than the formula of [3]. Instead 
of cutting along a system of edges of a graph embedded into a Riemann surface by a separat- 
ing curve, we consider a separating curve containing vertices only of the graph. To obtain a 
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cutting (gluing) formula for partition functions, we must perform a certain graph-theoretical 
constrution on the graphs resulting from the cut (introducing new subgraphs which we call 
core in the case of the dimer and target in the case of the Ising model). Our main gluing 
theorems are Theorems [TD] and of Section 121 

To work our way closer toward a limit formula, we need to understand how gluing works 
mathematically in the fermion conformal field theory. This is done in Section |H Even 
though there are both mathematical and physical approaches (e.g. [311 1221 [9]), as far as we 
know, there is no treatment in the literature which would include a mathematically precise 
definition of the data at a closed Riemann surface in terms of the Dirac operator, and the 
related cases of gluing. We write the Dirac operator on a Riemann surface in a conformally 
invariant form, and write down a definition and gluing properties of its determinant with 
respect to gluing along curves with antiperiodic spin structure. (The determinant is the 
square of the Pfaffian: we restrict to the case of the determinant because it is much simpler 
technically.) 

A rigorous discrete-to-continuous limit formula for one of the chiral partition functions of 
the chiral fermion on a torus is the topic of Section El Previous work on this was done by 
Cohn, Kenyon and Propp who work out a formula in the case of rectangular dominoes [6] 
(in fact, in the case of a parallelogram, it was worked out already by Kasteleyn [20], see also 
Ferdinand [10]). A heuristic approach in a more general situation is done in [8]. A general 
obstacle to deal with, which we haven't mentioned above, is the conformal anomaly, namely 
the fact that the elements U{T,) are defined only up to scalar multiple. (In the chiral sector, 
this propagates into the elaborate structure of a modular functor, some of which is described 
in Section O) To deal with this difficulty, one encounters a phenomenon in nearly all discrete- 
to-continuous limits in physics: regularization. Typically, one encounters a limit which does 
not exist, and is "regularized" by replacing it with a somewhat different but compellingly 
similar expression. In the present setting, we prove for arbitrary periodic critical embedding 
in a torus that the limit does not exist, and can be regularized by replacing, at one point, a 
set of parameters by their absolute values, after which the CFT limit is obtained. 

There is yet another complication in a rigorous approach to discrete-to-continuous limits 
in the case of the Ising model, which does not arise in the case of the dimer model: while in 
the case of the dimer model, the limit theory is the 2-dimensional fermion CFT, for the Ising 
model, the limit theory should be the 1-dimensional fermion [29], which, despite similar 
terminology, is considerably more subtle, and involves a construction of P. Deligne using 
super-central simple algebras over C. Whether natural discrete versions of these concepts 
exist is an open problem, which we discuss at the end of Section |6l 

2. Discrete gluing: the set-up 

Let G = (y{G), E{G)) be a finite unoriented graph consisting of vertices V{G) and edges 
E(G), in which loop-edges and multiple edges are allowed. We say that the subset E' C E{G) 
is even if the graph {y{G),E') has even degree (icluding degree zero) at each vertex. We 
say that M C E{G) is a perfect matching or dimer arrangement if the graph {y{G), M) has 
degree one at each vertex. Let £{G) denote the set of all even subgraphs of G, and let V{G) 
denote the set of all perfect matchings (dimers) of G. 
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We assume that an indeterminate Xe is associated with each edge e G E{G), and Xe may 
be evaluated in complex numbers. We define the generating functions for even subsets of 
E{G) resp. for perfect matchings, £g resp. Vg, as elements of Z[(xe)eeE(G)]: 

E'&£{G) eeE' 

'Pg{x)= Y1 U"^- 

Polynomial Vg{x) is also known as the dimer partition function of graph G. Knowing the 
polynomial £g is equivalent to knowing the partition function Z^"^ of the Ising model on 
the graph G, defined by 

where the Je (e G E{G)) are weights (coupling constants) associated with the edges of the 
graph G, the parameter /3 is the inverse temperature, and 

E n 

a:V{G)^{l-l} e={u,v}€E{G) 



The theorem of van der Waerden [33j states that Z^™^(a;) is equivalent to Sg{x) up to a 
change of variables and multiplication by a constant factor: 




Let us assume the vertices of G are numbered from 1 to n. If D is an orientation of G, 
we denote by A{G, D) the skew-symmetric adjacency matrix of D defined as follows: The 
diagonal entries of A{G, D) are zero, and the off-diagonal entries are 

A{G,D),^ = J2±^e, 

where the sum is over all edges e connecting vertices i and j, and the sign in front of Xe is 
1 if e is oriented from i to j in the orientation D, and —1 otherwise. As is well-known, the 
Pfaffian of this matrix counts perfect matchings of G with signs: 

Pfaf A{G, D)= J2 sign(M, D) JJ , 

M&V{G) eGM 

where sign{M,D) = ±1. We use this as the definition of the sign of a perfect matching 
M with respect to an orientation D. It is well known that we have for any pair of perfect 
matchings M, N 

sign(M, D) ■ sign(A^, D) = (-1) °^ ^"^^^'^ ^y^^"' 

An even- length cycle is even if it has an even number of directed edges in agreement 
with one way of traversal. 

We denote the polynomial Pfaf A{G, D) G Z[(xe)ee£;{G)] by Fd{x) and call it the Pfaffian 
associated to the orientation D. The following result is well-known: 
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Theorem 1 (Kasteleyn [20], Galluccio-Loebl [15], Tesler [32], Cimasoni-Reshetikhin [1]). 

// G embeds into an orientable surface of genus g, then there exist # orientations Di (i = 

1. . . . , A^) of G such that the perfect matching polynomial Vg{x) can be expressed as a linear 
combination of the Pfaffian polynomials F^- (x) . 

We call this expression the Arf-invariant formula, as it is based on a property of Arf 
invariant for quadratic forms in characteristic two. As far as we know, the relationship with 
Arf invariant was first observed in [1|. 

Let the graph G be embedded in a closed Riemann surface X of genus g. We denote 
H := F2) the first homology group of X with coefficients in the field F2. Recall that 

H carries a non-degenerate skew-symmetric bilinear form called (mod 2) intersection form 
and denoted by Quadratic form on {H, ■) associated to '■' is a function q : H ^ ¥2 
fulfilling 

q[x + y) = q{x) + q{y) + x-y 

for all x, y E H . 

We denote Q the set of quadratic forms on H over F2, its cardinality being 4^. 
Each quadratic form q E Q determines the signed generating function of even subsets of 
G by the formula 

(2) soA^) = i-^y^^"'^^ n 

E'e£{G) e&E' 

The generating function Sg{x) of the even subsets of edges admits analogous Arf invari- 
ant formula, where Sg{x) is expressed as a linear combination of its signed modifications 
£G,q{x),q e Q, see [25]. 

Analogously, partition function of the free fermion conformal field theory in (spacetime) 
dimension 2 on a closed Riemann surface X of genus g can be expressed as a linear com- 
bination of 2^^ functional Pfaffians of the Dirac operators, each of them corresponding to a 
Spin-structure on X. From a physical perspective, this was noted in e.g. [T], and rigorous 
mathematical approaches now also exist (see Sections IH E] below). 

Quite recently an extensive effort was invested in the understanding of the Dimer and Ising 
models, their criticality resp. asymptotic behaviour in terms of the conformal invariance 
and quantum field theory (conformal field theory in the critical case) of the free fermion and 
related theories (see e.g. [1], [7], [31], [5], [25], [8]). We continue in this effort and show 
that the dimer and the Ising partition functions share with the partition function of the free 
fermion a basic gluing property of their signed modifications. 

We present the gluing formulas for finite graphs embedded in double torus = T^T, 
regarded as a connected sum of two 2-tori (elliptic curves) T. The restriction is done for the 
sake of simplicity, i.e. analogous formulas can be proven for closed Riemann surfaces of any 
genus. 

2.1. Exterior algebra of gluing variables. Let us consider a graph G embedded in a 
double torus T^. We denote by C a simple closed curve (a cycle) on with the following 
properties satisfied: 

(1) \ C has two connected components, each of which is a torus with a disc bounded 
by G cut out. 
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(2) C intersects embedded graph G in an even number of vertices and in no inner point 
of an edge. 

(3) G has an even number of vertices in each component of \ G . 




Figure 1. The double torus 

Let v^, . . . ,v'', k even, be the vertices of G on C ordered along G. We consider the exterior 
algebra of 2k gluing variables 

e\ < el < ... < e\< e2 

ordered according to the orientation of G and the choice of a vertex on G. Then the vector 
space on the basis B consisting of all monomials 



(3) 
with 



e\\ < < ... < e ■ 



is naturally identified with the exterior algebra 

^25 •••5 ^2)- 

The dual exterior algebra 

^(^1, 62, 6]^, 62) = A(e2 ,6]^ ,...,62 ,6]^ ) 

is naturally identified with the vector space on the basis B* consisting of all monomials 



(4) 



It* 11* 



Let us consider Y G A(eJ,e2, ...,6^,62) written as a linear combination of basis elements 
Y = YlibeBV^^^ ^^"^ similarly Z G K{e\,e\, ...,6^,62)* written as Z = J2b'eB* ^b'b'. Then the 
scalar product < Y, Z > is defined by 

<Y,Z >=^ < ybb,Zb'b' >, 



where < yijb^zyV > is written in the following way. Let b = e*i^e*^...e*^ and b' = e^ 



...e- 



Then < b,b' >= unless m = i and ep = ap,ip = jp for each p, in which case < ybb, zyb' >- 
VbZb'- 

We first illustrate the gluing by writing down a gluing formula for the unsigned partition 
function £g{x)- This was done in a different way for the dimer partition function Vg{x) in 



2.2. Discrete unsigned gluing. We denote by Gi,G2 the subgraphs of G supported on 
connected components of \ C along with the vertices of G supported on G. Hence, if v 
is a vertex of G embedded in G then the set of the edges incident with v is partitioned 
into two sets Bi{v), B2{v), the first belonging to Gi and the second to G2- 

We construct an auxiliary graph, called the target, as follows. For each i = l,...,k, 
let us introduce the path L{i) of length 2 containing the vertices v{,V2jVI and the edges 
{^2^1}- The target S is obtained as the union of paths L{i),i = 1, . . . , by identifying 
all k vertices vl,i = 1, . . . , k into one single vertex (see Figure [2] for the target S). The target 
S allows to construct the graph Gj from Gj and S by identifying, for each i = 1, . . . , k the 
vertices and v\. 

The edge- variables of S are the gluing variables. We denote by = e\, = 63 the 

gluing variables in Gf and by y2vi^vi^ = ^2*5 y2vi^v\ = ^^e gluing variables in Gf. 



The vector yi of variables of Gf consists of {xe)e£E{G) and the gluing variables {e\, el, e^, 62). 
The vector y2 of variables of Gf consists of {xe)e£E{G) and the gluing variables (cg*, e^*, ...,el* ,el*). 
The following theorem simply follows from the definition of scalar product in the exterior 
algebra. 

Theorem 2. 



Many results in this section exploit the structure of (7-graph, extensively used and explained 
in detail in, e.g. [24j, [25j and references therein. 

Definition 3. The highway surface Sg consists of a base polygon Rq and collection of bridges 
Ri, . . . , R2g, where 

(1) -Ro is a convex Ag—gon with vertices ai, . . . , a^g numbered clockwise, 

(2) Each R2i-i is a rectangle with vertices x{i,l), . . . ,x{i,4) numbered clockwise, and 
glued to Rq in a way that its edge [x{i, l),x{i, 2)] is identified with the edge [ci4(t_i)+i, a4(i-i)+2] 
and the edge [x{i,3),x{i,4:)] is identified with the edge [a4(j_i)+3, a4(j_i)+4], 

(3) Each i?2i is a rectangle with vertices y{i, 1), . . . , y{i, 4) numbered clockwise. It is glued 
to Rq so that its edge [y{i,l),y{i,2)] is identified with the edge [a4(j„i)+2, a4(j-i)+3] 
and the edge [y{i,3),y{i,4:)] is identified with the edge [a4(j_i)_|_4, a4(.j_i)_|_5] . Notice 
that indices are considered modulo 4g. 

There is an orientation-preserving map $ of E>g into the plane M?, which is an immersion 
up to the images of the bridges i?2i and R2i-i intersecting in a square for each i = 1, . . . g. 
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Figure 2. The target 5* 



Sg{x) =< £as{yi),£cs{y2) > . 



3. Discrete signed gluing 



Let us assume the graph G is embedded into a closed orientable surface X of genus g. The 
combinatorial model of X is §g unified with an additional disk i?oo glued to the boundary 
of Sg. By an isotopy property of the embedding, we may assume that G does not meet the 
disk i?oo and moreover, all vertices of G lie in the interior of Rq. We may also assume that 
the intersection of G with any of the rectangular bridges Ri consists of disjoint straight lines 
connecting the two sides of Ri glued to the base polygon Rq. This assumption is not really 
necessary, but makes all considerations below more transparent. The composition of the 
embedding of G into with the immersion $ yields a drawing of G in the plane M^. A 
planar drawing of G obtained in this way will be called special. Observe that double points 
of a special drawing can only come from the intersection of the images of bridges under the 
immersion $ of Sg into the plane. Thus every double point of a special drawing lies in one 
of the squares $(-R2j) n <l>(i?2i-i)- 

In what follows we employ this machinary to produce gluing formulas. 

3.1. The Ising model. As for the notation and basic results, we follow [25] . 

Definition 4. Let G he a graph embedded in and let e he an edge of G. By definition, 
the embedding of e intersects each bridge Ri in a collection of disjoint straight lines, whose 
number is denoted rj(e). For a subset A of edges of G we denote r{A) the vector of length 
2g, whose i-th component is r{A)i = J2eeA^ii^) ■ 

We observe that two even subsets of edges A, B belong to the same homology class in 
H if and only if r{A) = r{B) (mod 2). We pick a basis ai, 6i, . . . , a^, bg of H, where a^, bi 
correspond to the even subsets Ai, Bi satisfying r(Aj)2i_i = r{Bi)2i = 1 and all remaining 
components of r{Ai) and r{Bi) are zero. Note that each quadratic form in Q is uniquelly 
determined by its values on any basis of the underlying vector space. 

The following Theorem is a basic result of ( p5 j ) . 

Theorem 5. Let G be a graph embedded in a closed Riemann surface X of genus g and let 
q be a quadratic form on the homology group H. Let E' be an even subset of edges of G and 
let r = r{E'). Then 



The importance of Theorem [5] is related to the following observation. Let us consider 
an embedding of G in the highway surface §2, for simplicity considered of genus g = 2. 
A quadratic form q E Q naturally determines two quadratic forms gi,g2 on the torus by 
= 1 (liipi) = qipi) resp. ^2(01) = ^'(02), 92(^1) = ^(^'2) together with trivial values 
on the other basis elements not written explicitly. 

The signed version of Theorem [2] is of the following form. 

Theorem 6. We have 




where 
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1=1 



(5) 



S{G,q,x) =< S{Gf,qi,yi),S{Gl,q2,y2) > ■ 



Proof: Using Theorem [2] the proof is easily reduced to 
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Claim. Let E' be an even subset of edges of G. For i = 1,2, there is exactly one extension 
of E' n Gi to an even subset of Gf by edges in the target. Moreover, for arbitrary q E Q, 

^_iy{lE']) ^ ^_iyi{lEi]) i^_iy2{lEi,]) _ 

The Claim follows directly from Theorem O and definitions of qi,q2- This finishes the proof 
of Theorem [61 

□ 

3.2. The dimer model. In this subsection we present a gluing formula for the signed 
dimer partition function of a graph G embedded in the closed Riemann surface of genus 
2. Analogously to the case of the Ising partition function, we first introduce an auxiliary 
graph called the core and denoted 5". The construction of the core follows the procedure of 
reduction of the Ising model to the dimer model, see |13j . 

We recall that in Subsection 12.21 we denoted by Gi,G2 the subgraphs of G supported on 
connected components of \ C together with the vertices of G supported on C. If f is a 
vertex of G embedded in C then the set By of the edges incident with v is partitioned into two 
sets Bi{v), B2{y), the first belonging to Gi and the second to G2- For each i = 1, . . . ,k, L{i) 
denotes the path of length 2 consisting of the vertices v\,v\, v\ and the edges {^3^2}, {^2^1}- 

To construct the core S' , we first consider the path consisting of Qk new vertices Ui, . . . , -Ugfc 
and edges {ui, Wi+i}, i = 1, . . .Qk — 1. Then we add the edges {u^j-21 U3j},j = 1, 2k and 
identify, for i = 1, . . . ,k, the vertex 'U6i-4 with (see Figure Elfor an orientation of the core). 

Finally, for j = 1, 2 we construct the graphs Gj from the union of Gj and S' by identifying, 
for each i = 1, . . . , k the vertices f * and v{. 

The relevance of the core stems from the following observation, see e.g. |13] . 

Observation 7. Let A C {v[, i = 1, . . . ,k}. Then S' has a matching completely covering all 
vertices of S'\A if and only if \A\ is even. Moreover, if\A\ is even, the matching completely 
covering S' \A is uniquelly determined. 

The edge-variables of the paths L{i) are, as in the construction of the target, the gluing 
variables. For i = 1, . . . , fc, we denote by Zi^i^yi^ = e\, z^^^^^^^ = 63 the gluing variables in Gf 
resp. 2:2^1^. = e'2*, ^2«^^j = el* the gluing variables in Gf'. 

We let vector zi of variables of Gf consist of {xe)eeE(G)y the gluing variables (e}, el, e\, 62) 
and the remaining edge variables be set to 1. We let vector Z2 of the variables of Gf' consist 
of 

{xe)e&E{G)^ the gluiug Variables {62*, ef*, ...,el*,el*), and the remaining edge variables be 
set to 1. 

We are now ready to discuss gluing formula for unsigned partition function. The follow- 
ing Theorem simply follows from the definition of the trace in the exterior algebra, using 
Observation [71 

Theorem 8. 

Vg{x) =<V^s.{z,),V^s.{z2)> 

In the next step we consider the signed case. We first modify G by adding the edges 
{t>*,f*+^},i = 1, . . . ,k — 1 and let the variables associated to these edges be equal to 0. 
Hence, adding this set of edges does not change the signed dimer partition function of G. 
Let us denote the modified graph by G'. 
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Figure 3. The orientation of the core S' 

A result of [32] (see also recent [21]) can be summarized, in the case of the double torus, 
as follows: 

Theorem 9. Let G be a graph embedded in a highway surface §2- Let D be an orientation 
of G so that for each bridge R, the immersion of Rq U R into the plane has the property 
that each inner face has an odd number of edges oriented clokwise in D. Let M be a perfect 
matching ofG. Then 

sign(M,D) = 

where c(M) denotes the number of self-intersections among edges of M inside the special 
drawing of G in the plane. 

Let G be the graph embedded in the double torus T^. We consider the embedding of G' 
in the highway surface S2 determined by an embedding of G in the closed Riemann surface 
of genus 2 . Note that the new edges attached among the vertices of G form a chord of 
the base polygon Rq, which separates the two pairs of bridges. We recall that f^, . . . is 
the order of the vertices of G on C along this chord, and that k is even. 

It is not difficult to observe that there is an orientation Dq of G' fulfilling the properties 
of Theorem [9], where the orientations of the edges {t;*, v"^^^}, i = l,...,A; — 1, form a directed 
path of odd length from to v^. 

Let D be an orientation of G. We denote by Di the restriction of D to Gi and by D2 the 
restriction of D io G2- We consider the orientation of the core S", indicated on Figure [31 By 
slight abuse of notation we denote it also by S' . If D is an orientation of G then we denote 
by Df , Df the orientations of Gf , Gf induced by D and S' . 

Theorem 10. Let D be an orientation of G. Then 

Fd{x) =< F^s.{z{),F^s,{z2) > . 

Proof: Let M be a perfect matching of G. It follows from Observation [7] that for i = 1,2, 
there is exactly one extension Mj of M fl Gj to perfect matching of Gf by the core edges. 
Theorem [TU] follows from Theorem [H] and 
Claim. 

sign(M, D) = sign(Mi, Df) ■ sign(M2, D^). 

We first show that the Claim holds for the orientation Dq introduced above. Since M is a 
subset of edges of G, we consider Dq restricted to G only. We have by Theorem [9] 

sign(M,Do) = (-l)'^*'\ 
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where c(M) denotes the number of self- intersections among the edges of M in the special 
drawing of G' in the plane. From the construction of G' in Subsection 13.21 we have 

c(M) = ci(M) + C2(M), 

where Ci(M) denotes the number of self- intersections among the edges of M that belong to 
the first pair of the bridges of the highway surface §2, and C2(M) denotes the number of 
self-intersections among the edges of M that belong to the second pair of the bridges of the 
highway surface §2- 

Finally, the orientations Df , are constructed from Dq in such a way that they also 
satisfy the properties required by Theorem [91 Hence 

sign(Mi,Df) = (-ir(*^i), 

(6) sign(M2,Df) = 

Summarising, the Claim holds for the orientation Dq. An arbitrary orientation D of G is 
obtained from Dq by reversing orientation in a subset S{D) of the directed edges of Z^oi so 

sign(M,D) = sign(M,Do)(-l)'^^^^''^' 

and 

sign(Mi, Z^f) = sign(Mi,Df 

(7) sign(M2,Df ) = sign(M2,Df 

This finishes the proof of the last Claim and at the same time of Theorem [TOl 

□ 

4. Fermions on compact surfaces I 

In this section, we present some facts about the free fermion conformal field theory in the 
context of the desired limit from the discrete to the continuous case. From the physics point of 
view, this material is more than two decades old, see e.g. [H |3l]. A mathematical approach 
was outlined in [31], and further developed in [22| [18] . Nevertheless, the mathematical 
literature as it stands is deficient in an important point (explained below), which happens 
to be the key to addressing the limit questions we are interested in. Treating this gap is the 
main purpose of the present section. 

The point is that in the mathematical treatments of CFT, dosed Riemann surfaces (or 
compact Riemann surfaces without boundary) play a special role, and a technical complica- 
tion. It is a common technique, in fact, to consider only (compact) Riemann surfaces without 
closed connected components, and argue that data at closed surfaces can be recovered by 
cutting out a holomorphic disk, and giving compatible isomorphisms of the data as the disk 
varies. This is the approach taken in pT| [22] . 

For our purposes, however, this is unsatisfactory. The reason is that from the physical 
point of view, the case of a closed surface is in fact more fundamental ([34j), because it is 
where contact with quantum field theory is made. The data we are supposed to have in 
this case is the determinant (resp. Pfaffian) of the Dirac operator, which is analogous to the 
discrete data we studied in sections [H [H 
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The purpose of this section, then, is to explain how to define the determinant precisely, and 
to describe gluing from the data on a surface with boundary as described in [311 [22] to the 
Dirac operator data in the closed surface case. For simplicity, we treat only the determinant, 
which corresponds to the free fermion of spacetime dimension 2, and not the Pfaffian, which 
is more complicated. Also for reasons of simplicity, we discuss only antiperiodic boundary 
conditions. The role of all spin structures will be reviewed in more detail in Section [21 
below. A discrete analogue of the 'boundary to no boundary' gluing case will be, in fact, 
exhibited in Section [5] below, where we give one prototype of a rigorous regularized CFT 
limit computation. 

Let E be a compact Riemann surface. A spin structure on S can be identified with a 
choice of a holomorphic line bundle il^^'^ on S and an isomorphism 

(We identify a holomorphic bundle with its sheaf of sections M; we denote by M the space 
of global holomorphic sections; M^^, M^^ denotes the corresponding notions in the smooth 
category.) The bundles of spinors on S are A~ = and = ^l'^^^'^ . 

The Dirac operator 

ordinarily depends on metric (i.e. is only defined on a spin Riemann manifold), which is 
undesirable for the purposes of CFT. However, a modified conformally invariant construction 
can be obtained as follows: We define a pairing 

(8) IP : nli^ ®c -> C 
in a local holomorphic coordinate z as follows: 

(9) fdz^^^ ® gdz^/^ ^ J ^gdzdz. 

To see that ([9|) is invariant under holomorphic coordinate change, let z = z{y). Using the y 
coordinate, we compute 

fdz^' ® gdz'/' = f l^^j dy'l^ ® g j dy'l' ^ 

f df dz f df dzdz f ^ ^- 

^ / T^9-rdydy= / -^9-r-p.dydy = / —gdzdz. 
J ay dy J oz dy dy J oz 

Also note that by Stokes' theorem, Ip is an antisymmetric pairing, since 

^9 + dzdz = d{fgdz). 

One defines (up to scalar multiple) the partition function of the 1-dimensional chiral fermion 
on E as 

The partition function of a fermion with both chiralities is 



(10) 



where Ip is the construction analogous to ij) with "holomorphic" replaced by "anti-holomorphic" 
and a{a) is the Arf invariant. 

To define the Pfaffian, one proceeds as follows. For an antisymmetric Fredholm operator 

where if is a Hilbert space (we are not in this situation, but may achieve it by suitable 
Hilbert-completion), there is a Pfaffian line Pf{^) (well defined up to canonical isomor- 
phism), and a Pfaffian element 

(11) p/(<|.) e P/(<l>) 

(see [30], Chapter 12). For simpicity, we will describe only the square of (ITT!) . 

(12) det{^) G Det{^) for ^ : Hi ^ H2 Fredholm of index 0. 
Choose an invertible operator J : Hi ^ H2 such that 

T = J- $ 

is trace class. (For any compact operator K : Hi ^ H2, there exist unique numbers si > S2 > 
... > 0, lims„ = and bases (e„), (/„) of Hi,H2, respectively, such that K{en) = Snfn- The 
numbers s„ are called singular values. K is trace class (resp. Hilbert- Schmidt) if ^ |s„| < 00 
(resp. Y^i^nf < 00).) 
Now one defines 

Det($) := ((J)) 
(=the line with free formal generator (J)), and 

det{^) := det{I - TT^) ■ (J) G Det{<^) 
(the determinant is defined as a number because TJ~^ is trace class). For another choice 

$ = / - T' 
with J' invertible and T' trace class, we have 

J'J-^ = I + {T' -T)J-\ 

so the canonical iso 
is 

{J') ^ (J) ■ det{I + {T - T)J-'). 

One has det{^) = unless $ is invertible. To define det{^) as a number, one must specify 
a generator of Det{^). This is called regularization. 

Let us turn back to Pfaffians. In our situation, this amounts to picking a section of 
the pullback of the bundle Pf{Ipj^) on the Teichmiiller space. This can be done using the 
structure of CFT. 

For a Riemann surface S with spin structure with analytically parametrized boundary 
such that the induced spin structure on each boundary component is antiperiodic, there 

1 /2 

is a non-degenerate quadratic form on the space Qqj^ of complex- valued 1/2-forms on the 
boundary given by 

B{uj,v) = / ^V- 

JdT. 
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For each boundary component, 1/2-forms which extend to the unit disk holomorphically 
form a maximal isotropic subspace of (Q^Ji). One defines 

(13) -Hsi = iA{nf rr 

(here ?^ denotes Hilbert completion,) and 

Has := (g)(H5i)*®(g)H5i 

where ® denotes Hilbert tensor product, and the tensor products on the right hand side are 
over boundary components with inbound resp. outbound orientation (the orientation of the 
boundary component of the standard disk D parametrized by the identity is outbound). 

In fact, more explanation is in order: In the present situation, note that the symmetric 
bilinear form B together with the Real structure (complex conjugation) specifies an inner 

product on Q^J^ (hence on flg^^) which is the specific one we use in forming the Hilbert 
completion, (and hence also defining, say, the inner product on A((fi^(^)~) etc. 
Now the inclusion 

(14) (fif )- c (n'J^r 

where the second map is orthogonal projection is an inclusion of a maximal isotropic subspace 

1 /2 

which is restrictedin the sense that its orthogonal projection to {{^q-^ )~)^ is Fredholm, while 

1 /2 

the orthogonal projection to the orthogonal complement ((^^gs)^)^ Hilbert-Schmidt. 
In this situation, fl]^'^ specifies a line Pf{fl]^'^) and a linear map 

(15) Pfi^T) -> ^OE. 

Now when S is a closed Riemann surface with spin structure, define 

(16) Pf{nf):=Pf{l/)^). 

Lemma 11. Let 11 be a Riemann surface obtained from S by gluing some inbound boundary 
components to outbound boundary components in a way so as to preserve parametrization. 
Then this data specifies an isomorphism 

(17) P/(fif)^P/(fif) 

and a commutative diagram where the rows are given by CFT structure 

PfinU') -Ha^ 



pfiSif) — 



tr 



as- 

(In the case when S is closed, the bottom right corner is C, and the bottom row is defined to 
be pf{Ip)-\) 

We will show only the construction of the map ( ITTll in the case when S is closed (i.e. 
where the definitions of the two Pfaffian lines are different). Also, for simplicity, we will only 
treat the square, i.e. the analogous construction for determinant lines: the construction for 
Pfaffians is similar but more elaborate. 
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Denote by (fiQ'^^)^ the subspace of the space of smooth sections u of (il]m)s with the 
property that for any two boundary components in (9S which are to be glued in £, 



du:\ 



Cl 



duj\ 



C2- 



Denote the image of in E by 9S/2. Then we have a commutative diagram 



1/2 



'(9S/2 



(19) 



{{Qsrn )g) 



1] 



1/2 

aE/2 



where all rows and columns are short exact with the exception of the top row and left column, 
which are Fredholm (to reduce clutter, we omit the symbols for Hilbert completions). The 
top row differs from (IT^ by a trace class operator (note: ultimately, this difference gives rise 
to the eta-factor!), so it suffices to prove that for any diagram of Hilbert spaces and bounded 
operators 



(20) 



C 



A' 



D 



W 



A 



G 



B 



B 



where all the rows and columns are short exact except F and G which are Fredholm, we 
have a canonical isomorphism 

(21) Det{F) = Det{G). 

But note that in ( l20l) . we may write B = W/C, A = W/D, so both sides are obviously 
canonically isomorphic to 

(22) Det{i®j:C®D^W). 

The commutativity of the square (fT8l) is proved in [22]. □ 

Now for the case = Ag (the standard annulus in C with boundary components parametrized 
by z, qz and antiperiodic spin structure), the Lemma specifies a canonical isomorphism 

which specifies a canonical element 

(23) e Pfi^X). 

Considering the elliptic curve 11 = Eg over C, the image 



ig := iiig) e Pf{n%^ 
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1/2^ 



is the desired regularization. With respect to this regularization, the partition function is 
one of the two Jacobi theta functions involving odd powers of g^^^ (depending on the spin 
of the gluing), times the factor 



(24) n (1 - Q") = V{r)q^^^\ q = e'"'". 

n>0 

The remaining Arf invariant theta function is not accessible in this way, but may be 
obtained by cutting around a different simple curve in Eg of anti-periodic spin structure; 
after multiplying by the factor it is related to the remaining Arf invariant partition 

functions by a modular transformation. The Arf invariant 1 partition function is 0. 

In the next section, we will see that the partition function of the chiral fermion is, in a 
special case where the boundary can be easily modelled discretely, a "regularized limit" of 
partition functions of dimer models on critically embedded graphs, using a discrete analogue 
of Lemma [TT] 

In Section [2] Theorem [8|, we exhibited a discrete analogue of Lemma [11] for arbitrary (not 
necessarily critically embedded) finite graphs, embedded into a Riemann surface of genus 2. 



5. A LIMIT FORMULA FOR THE DISCRETE DIMER MODEL ON A TORUS 



We use here the setup of Kenyon [21]. Consider a bipartite planar graph with vertices 
colored black and white. The set of black (resp. white) vertices is denoted by B resp. W. 
A critical embedding into a plane is such that the vertices of each face lie on a circle, and 
all these circles have the same radius. One defines a square matrix d with rows and columns 
indexed hj B UW as follows: If Vi,V2 are not adjacent, then d{vi,V2) = 0. If w and b are 
adjacent vertices, w being white and b being black, then 

(25) d{w,b) = -d{b,w) 

is the complex number of length given by twice the distance of the center of the circle 
circumscribed around (either) face containing b, w to the mid-point between b and w, and 
direction pointing from w to b. This is the total Kasteleyn matrix, which we denote by 

K 
-K 

where the top (resp. bottom) set of rows corresponds to white (resp. black) vertices. We 
refer to K as the (discrete) chiral Kasteleyn matrix. (In ( l25l) . our convention differs from [21] ; 
where the matrix is symmetrical; however, antisymmetric matrices fit better with Pfaffians.) 

In view of the discrete field theory of Cimasoni-Reshetikhin, the material of Section H] 
above suggests a formula of the following form: Let F be a graph critically embedded into 
C/Z{l,r} such that [0,1] is antiperiodic. Let 6 be a black vertex and let fb{z) be the 
corresponding discrete exponential. Then one should have for the chiral Kasteleyn matrix 



(26) Det{K)^ n {l±h+r{z)) 

z-fb+i{z)=-fb{z) 
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where the sign is the spin structure on [0, r]. 

This formula is further motivated by the resuh of Mercat [26] that exponentials generate 
additively the vector space of discrete holomorphic functions; rewritten in the basis of discrete 
exponentials, the chiral Kasteleyn matrix becomes diagonal, so its determinant becomes the 
product of the diagonal terms, analogously to formula fl26l) in Section H] above. There are 
however two caveats to applying Mercat 's theorem: first of all, Mercat is considering a 
slightly different setup, not assuming that the graph F is bipartite. He works with tilings by 
rhomboids; the graph made by the edges of the rhomboids is naturally bipartite, with the 
vertices of the graphs F, F* being the vertices of the two different colors. Second, Mercat 's 
result applies to tilings of simply connected subsets of C, and he does not specify how 
many discrete exponentials are needed to get a basis of the space of discrete holomorphic 
functions; in order for the formula fl26l) to be strictly correct, we would have to prove that 
discrete exponentials f'?{z) for such z that 

(27) = -Mz) 

precisely form a basis of discrete holomorphic functions on F tiling Ag twisted by the antiperi- 
odic spin structure on [0, 1] , g = 27rir. One precise statement (for the case of a rhomboid 
graph) is as follows: 

Proposition 12. Consider a bipartite graph F critically embedded into the torus Ej. = 
C/(l,r) (Im{T) > 0) in such a way that there exist complex numbers ai,...,a2n, &i,---,&2m 
of equal absolute value such that ai + ... + = 1, 6i + ... + 62m = and the corresponding 
rhomboid graph on F U F* has edges 

(ai + ... + Ofc + 61 + ... + hi, ai + ... + flfc+i + bi + ... + be), 

(ai + ... + Ofc + 61 + ... + be, ai + ... + Ofc + 61 + ... + 6^+1) 

(here we consider the subscripts as elements ofZ/m, Z/n). Assume further that the vertex 
is black, and assume that the spin structure of Er around [0, 1] is antiperiodic. Then the 
Kasteleyn matrix of F satisfies 

(28) pf{Kr) = 2™ fi (1 T fr{z,)) 

where fb{z) are Kenyon exponentials [21], and Zj ranges over all complex numbers for which 

(29) /i(z,) = -1, 
and the sign depends on the spin structure of [0,r]. 

Remark: This is not as easy to generalize as one may think. In particular, in case of 
more general Riemann surfaces, a straightforward generalization of this formula is false. A 
heuristic explanation may be extracted from Section H] above: In general, we are unable 
to identify functions with 1/2-forms in boundary behavior, and transformation rules for 
1/2- forms would have to be modelled discretely. 
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Proof: Because of continuity, it suffices to prove the formula for generic values of Oj, bj. 
Consider for k = 0, ...,m — 1 the sets of black vertices of F 

Bk = {bi + ... + b2k + ai + ... + a + 2i\i = 0,...,n - 1}. 

and the sets of white vertices of F 

Wk = {6i + ... + b2k+i + ai + ...a2i+i \ i = 0,...,n- 1}. 

The Kenyon exponential function fb{z) with base point is on the black vertex 

b = ai + ... + + bi + ... + b2k 

equal to 



n 



n 



,=1 i-^&2. 1 



za2j 



Recall that we denoted by zi,...,Zn all numbers satisfying (!29|) . Consider functions 
where 

fb,k{zj) = fb{zj) when b E Bk, 
= when b ^ B^.. 

We claim that for generic values of Oj, bj, form a basis of the space of all functions 

5 — C. In effect, it is true for = ^ by direct computation, and hence it is true generically 
(since non- degeneracy occurs on a Zariski-open set). 
Now by direct computation. 



(30) 
where 



dfi^kizj) = g7,k{zj) - g7^k-i{zj) for k = 1, 
dfiflizj) = g?,o{zj) =F fT{zj)g?,n{zj) 



m — 1 



g,,k{zj) -.W 



are functions defined as follows 
(31) 



gbi+...+b2i+i+ai + ...+a2i+i,k\Zj) '.— 

/fei + ...+b2£+ai + ---+«2i+2,fc(%') ^ fbx+...+b2i+ax+...+a2i,k^Zj). 



Note that for each j, the matrix (1501) in the basis (/?, k{zj), g7^k{zj)) is 



V TfAzj) 1 
which has determinant 

^TfriZj). 

But now the function h7^k '. W ^ C which is defined by 

hbi + ...+b2e+i+ai + ...+a2i+i,k{Zj) := fbi+...+b2i+ai+...+a2i+2,k{Zj) 
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plays a symmetric role to f?^k which respect to exchanging d for d, and the base change 
between g and h is conjugate to the matrix 

/ 1 -1 \ 

1 -1 

1 -1 



V 



1 



1 



which is 2. 



□ 



To try to obtain a limit over tilings with decreasing length of rhomboid edge, it further 
seems natural to use a recipe of Mercat |26j how to take the limit in order to make discrete 
exponentials converge to ordinary continuous exponentials: subdivide each rhomboid into 4 
equal rhomboids with parallel edges, and repeat this procedure. Even in simple examples, 
however, we see that this procedure when applied to (126|) . will not produce a convergent 
limit in the naive sense: the limit has to be "regularized" (which is not surprising, given the 
usual physical context of lattice regularization) . For example, if the segment [0, 1] is tiled 
into equal segments by the edges of F, the equation (127|1 becomes 



26|) becomes a product of quantities of 



l + z/N 

which implies that z is on the imaginary line. Then 
the form 



(32) 



lie- 



so given Im{T) > 0, for Im(z) > 0, the second summand (l32ll will have absolute value > 1. 
The number of such terms will increase with , as solutions of ( 1271) in increasing bounded 
regions will approach {2k + l)7ri, G Z. Therefore, the expression f l2B]) in this case cannot 
converge as N ^ oo. The "limit" has to be taken in the sense of 



n (i-g"+^) = ± n (i-g-- 

neZ neN 



n (1-?^ 



nGN 



(the right hand side being the desired ratio of a theta function by rj), based on the "calcu- 
lation" 

1 3 5 



-1-1-5-2-5-... 

2 2 2 



:i + 2+3+...)~-lc(-i)+c(-i) = --^- 



In this special case, the regularization can be achieved by multiplying (!26|) by 



(33) 



where the product is over solutions of (1271) with Im{z) > 0. This regularization indeed 
produces the desired ratio of 9 and t] in the hmit, multiphed by the factor g^/^^. 

Therefore, in some sense, for a general graph, the question entails finding a regularization 
procedure which would generalize (133|) . Again, one precise statement can be obtained as 
follows. Define 

n 

■.pd{Kr):= n {l-fr{t\z,\)). 

Proposition 13. Suppose graphs Tj. are as in Proposition m, and that the rhomboid graph 
of Ffc+i is obtained from the rhomboid graph of by subdividing each rhomboid into 4 
congruent rhomboids with parallel edges. Then 



hm :pf{Kr,):= U {iTq- 

fc— >oo j=0 



Proof: We must study solutions zj of the equation 

n (l + 



We rewrite this as 



N 



(34) E (ln(l + ^)-lnfl-^))=^vr., fceZ. 



- . N J V N JJ N 
Put 

n 

fit) := E (ln(l + a2i-it) - ln(l - a2^t)). 

i=l 

Then 

/(O) = 

and / is analytic in a neighborhood of with non-zero derivative at 0. Hence, the same 
holds for the inverse g = f~^. From this, (substituting t = j^), we have the following 

Lemma 14. For every 6 > there exists a c> such that for every N, when 

(35) 1^1 < cN, 

then is in the domain of g, and all the solutions to [34\ ) satisfying l[3^) are of the form 

6k 

(36) Zk = i2k + l)Tci + tk, \ek\ < —. 



N 
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□ 



Now consider an inequality of the form 



(37) 



(j-l/M 



< 



l-q 



\z\/2-K 



1 _ q\k+in\ 

where C > 1 is some constant. ( 1371) holds when 



< 



g|fc+l/2||| < 



InC 



or 



\^^\z\/2-K _ ^|fc+l/2|| ^ 



InC 
2M' 



This follows whenever 
(38) 



\z\ 
2^ 



(|fc+2l) 



k 



for some constant > dependent on C, since g* has a bounded derivative in t > for 
< |g| < 1. 

Now set M := By the Lemma, there exist constants A, i? > 0, (we may assume 

A <1) such that ( l38l) (and hence (1371) ) holds for some solution 2;^ whenever 



(39) \k\ < AViV, 
and 

(40) \zj\> B^/N ioi a\\ ] e S 
where S is the set of all other indices j. Now note that 



(41) 



i65 



0{ E (|g''^'/'"|) = 0(Ar|g|^^/2-). 

ie5 



The limsup of the right hand side of ( 1411) with — )■ 00 is 0. Thus, the solutions ( HOl) may 
be neglected. Also, the quantity 



n (1 _ gl'=+V2|)2 
fc>AVJV 



approaches as — )■ 00, and by ([3 

1 



< 



n 



\k\<AVN 



1 _ gl^l/2^ 



1-g 



(|fc+l/2|) 



as M = yN and A < 1. Since C > 1 was an arbitrary constant, we are done. 



□ 
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6. FERMIONS on a RIEMANN surface II: A BRIEF REVIEW OF CONFORMAL FIELD 
THEORY, WITH SPECULATIONS ABOUT GRAPHS 

In physics, it is "known" that the critical Ising model on a lattice converges to the 1- 
dimensional free fermion theory [21], cf. [7], while the critical dimer model converges to 
the free fermion of space-time dimension 2 ([H])- To obtain rigorous mathematical theorems 
in these directions, we need to take into account the detailed structure of the fermionic 
CFT's, and model it in graph theory. Such mathematical treatments of conformal field 
theory now exist, (cf. [311 [22] from the Segal point of view, [9] from the vertex algebra 
point of view), but the structure involved is somewhat complicated. Because of this, we 
take time in this section to survery these structures in more detail. Some comments on 
their potential impact on the graph theory side of the story will be made at this section's 
conclusion. For simplicity, throughout this section, we will consider complex CFT structures, 
as real structures complicate things further. 

A big part of the story of fermion CFT's is bosonisation. For a bosonic CFT, we are 
supposed to have a real-analytic line bundle L on the moduli space of Riemann surfaces S 
with analytically parametrized boundary, a complex Hilbert space H and, writing Hg-£ for 

when S has m inbound and n outbound boundary components (we take C when S is closed), 
a map 

(42) ^ Hq^ 

satisfying gluing axioms analogous to Lemma [TT] (cf. [311 [181 flT])- 

A genus g Riemann surface Eg can be cut along g non-separating disjoint real-analytic 
simple curves Oi, into a genus surface A with 2g boundary components (let, for future 
reference, a "dual cutting" by real-analytic simple curves bi, ...,bg.). Gluing a standard disk 
on each of these boundary components, we obtain a sphere. Since the 2-sphere possesses a 
unique conformal structure, this specifies an element 

i'ai,...,ag e Lsg, 

hence, by f H2|) . a complex number. This assignment 

(43) Z{ai,...,ag) eC 

is called the genus g partition function. For g = 1, one can take the standard annulus Ag G C 
with boundary components and qS^ parametrized by the functions z, qz, respectively, 
which gives us the partition function Z[q). This function is a series in 

(44) 

(0 < ||g|| < 1) with positive integral coefficients where a — 6 G Z, and the coefficient of 
(jH]) counts the number of states of H of weight (a, b). In other words, it is the trace of the 
"grading operator" which, at weight {a,b), is given by multiplication by q°'q'^. 

There is a special class of CFT's called rational conformal field theory, abbr. RCFT. (One 
rigorous mathematical approach was reached by Huang and Lepowsky [19], although it 
doesn't use the language of Hilbert spaces, but the structure of vertex algebra. An approach 
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via Hilbert spaces was developed in [HI [TT]. To date, the two approaches have not been 
fully unified, and each has certain desirable features.) In this case, we have a finite set A 
of sectors with a distinguished element G A and an involution (?)* : A — )■ A, 0* = 0, and 
decompositions 

AeA 

Z{q) = Y,ZMZAq). 
AeA 

(Sometimes, elements of A are referred to as labels, an indexing set for sectors. We choose 
here to identify each label with its corresponding sector.) The precise explanation of (H5|) is 
that we have certain line bundles M\ on the moduli space of all annuli, and a positive-definite 
Hermitian pairing 

(46) 0M,®Mv ^L, 

AeA 

and we have 

Zx{q) : MxA, ^ C 

so makes sense. 

This is actually the tip of an iceberg called chiral conformal field theory: One is required 
to have a finite-dimensional complex vector space M^S for each Riemann surface ^{x^) with 
parametrized boundary components labeled by the Aj's, which are supposed to satisfy both 
fH6|) and gluing axioms [TH [11]. The most important part of the gluing structure is an 
isomorphism 

(47) ME(,,) = 0ME(,,,,,,) 

AeA 

where S is obtained from S by gluing an inbound and an outbound boundary component: 
the sum on the right hand side is over all possible labellings of the two additional boundary 
components by the same label. Additionally, there are to be linear maps 

(48) MS(,^) ^ i/aE,,^) 

satisfying appropriate gluing axioms ([HI [Hj). It is worth mentioning that for each A G 
A, the line bundle MxA specifies a holomorphic C^-central extension of the semigroup of 
annuli with analytically parametrized boundary components of opposite orientations; these 
central extensions do not depend on A, and are characterized by a single number called 
the central charge and denoted by c - see [31]. (Note: there is an infinite-dimensional 
space of annuli with analytically parametrized boundary components modulo the relation 
of conformal equivalence preserving the boundary parametrizations; if we restricted only to 
the standard annuli Aq, we would get the semigroup Co<5<i which does not have non-trivial 
central extensions.) 

Models of the fermion CFT which satisfy the above axioms are referred to as bosonised 
fermions. From the point of view of string theory, a fermion of spacetime dimension d has 
central charge d/2. The theory differs depending on whether d is odd or even. For d = 1, 
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the theory can be constructed, for example, as a coset model of the level 2 WZW-model for 
SU{2) ([7], 18.5.1, /c = 2), or as a non-supersymmetric minimal model [7], 7.4.2. 

Generally, for d odd, the d-dimensional bosonised fermion theory has three sectors 

A = {NS+, NS-, R}, A* = A, = NS+. 

(R stands for Ramond and NS stands for Neveu-Schwarz.) Its partition function is deter- 
mined by letting ZjV5'+(g), Z^s-io) be the integral (resp. integral +l/2)-dimensional terms 
of 

n>0 

and 

ZR{q) = 2^''^^\[{l + q-). 

n>0 

The power 2 in the formula for Z^^q) is significant, it expresses the fact that the bottom 
weight space of Hr is naturally the spinor module of the complex Clifford algebra on a space 
of dimension d. The dimension of the modular functor on a Riemann surface of genus g is 

(49) ^(4^ + 2^), 

which is equal to the number of spin structures of Arf invariant 0. 

The bosonised fermion of even spacetime dimension d has sector set 

A2 = {NS^,NS-,R+,R,}. 

We have 

R\ = R_. 

Here the labels we chose express the operations on the sector Hilbert spaces of the 1- 
dimensional fermion. In particular, Hr^Hr decomposes into two state spaces, according 
to the decomposition of the complex spinor on M?. Its partition function is determined by 
letting Zns+{.<1)i Z^s-il) be the integral (resp. integral +l/2)-dimensional terms of 

n(i + g"^^)', 

n>0 

and 

ZnM = 2''/'llil + qn. 

n>0 

The reason why there are now two Ramond sectors is that there are two non-isomorphic 
complex spinors in even dimension. The dimension of the modular functor on a Riemann 
surface of genus g is 

(50) 4^ 
which is equal to the number of all spin structures. 

For d = 2, the bosonised fermion theory turns out to be isomorphic to the conformal field 
theory associated with the lattice 2Z C M ( |T2l IT8]). This phenomenon is known as the 
boson-fermion correspondence. Expressed this way, the labels become 

{0,1/2,1,3/2}, (1/2)* = (3/2). 
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The chiral partition functions are 

(51) Z,{q) = J2q'^'^'^'^'l[{l-qn~\ 

We see that regardless of the parity of dimension, the partition function of the fermion of 
dimension d is 

\n>0 n>0 

(52) 

2'^+^(gg)'^/^6 J](l + g")(l+r) 

n>l 

Note that for = 2, this is indeed equal to 
(53) 

■ n (i-?")-^(i-r)-\ 

m,n>0 

which is what we obtain when we combine the chiral partition functions (!5T]) . Note that 
in ( |53|) and (!52|) . in fact the sums of the first two summands are equal by Jacobi's triple 
product identity, the last summands. 

How is it possible that the dimensions of the modular functor f H^ and f lSU]) in the even and 
odd- dimensional case are different, even though the partition function ( 152|) has a uniform 
expression? The answer is that in the even-dimensional case, when the spin structure along 
some of the curves is anti-periodic, the two Ramond sectors arising will have the same 
partition function, and hence their sum can be grouped in the power of 2 coefficient. 

The genuinely fermionic approach to fermion RCFT uses another, more complicated, 
axiomatization. This time, the chiral parts are chiral super conformal field theories. (This 
is a different concept from a supersymmetric conformal field theory, where both a bosonic 
and a fermionic part are present.) Chiral super conformal field theories are axiomatized 
similarly as chiral conformal field theories, but the Riemann surfaces S involved come with 
spin structure. We then have a decomposition 

A = Ajv^nA^j, 

expressing which labels apply to antiperiodic and which to periodic boundary components. 
Additionally, we deal with super-modular functors, which means that each MS(Aj) is a Z/2- 
graded vector space, as is each sector state space 

The maps fH5]) are required to be graded, and the usual sign convention 

(-1)^^- 

is applied with interchanging elements of degree i, j. We now have chiral partition functions 
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where the superscript expresses the spin structure in the radial direction of the annulus. 
When factoring a CFT to fermionic chiral factors, the analog of formula ( H5|) at both chiral- 
ities becomes 

(54) Z{q) = \ Zr{q)Z^,,'(s) + E (^)^f*®) • 

\AeA AeA / 

For a Riemann surface of arbitrary genus g where the summation is over spin structures on 
the Riemann surface E, and over a product of the sets of periodic (resp. antiperiodic) labels 
matching to the spin structure of each of the curves a,, 6,, and the coefficient before the 
sum is 1/2^. (This can be thought of as an averaging over spin structures, or, in analogy 
with string theory, a "GSO projection".) When there is exactly one periodic and one anti- 
periodic label, one can show that the super-modular functor is always one- dimensional for 
each spin structure on a Riemann surface, and is invertible in the sense that its tensor 
product with another super-modular functor is the trivial (unit) super-modular functor. In 
this case, the summation is precisely over the set of spin structures on S. This, in fact, is 
the CFT explanation of the decomposition of the discrete partition functions of the Ising 
and dimer models into summands corresponding to spin structures, with signs given by the 
Arf invariant (see e.g. formula [2] above); those are the summands whose gluing we study in 
Sections [2] and [31 There is one substantial difference in the "CFT limit" , however: one can 
show that the summands corresponding to spin structures of Arf invariant 1 are always equal 
to 0. From a physics point of view, this can be justified by noting that the contribution of 
such terms would come with a negative sign, which is physically impossible. A mathematical 
argument can also be obtained along the lines of [22] . 

It is also possible to talk about a fermionic conformal field theory with both chiralities, 
which consists of two super-Hilbert spaces and H^s (-R corresponds to the periodic 
and NS to the anti-periodic spin structure on S^), a real- analytic line bundle L on the 
moduli space of Riemann surfaces with spin structure and with (analytically) parametrized 
boundary components. The conformal field theory then specifies, again, a map 

Ly; — ^ Hqy 

where Hgj] is, again, a Hilbert tensor product of copies of Hr, H^s and their duals, depending 
on orientation and spin structure of each boundary components. The gluing axioms are 
analogous to the bosonic case (Note: no super-structure is present on L, or in other words, 
it is assumed to be "even".) 

Given a fermionic chiral RCFT, the spaces Hji and H^s are given by 

0i^A®^A* 
A 

where the sum is taken over R resp. NS labels. For a closed Riemann surface S, one then has 
one partition function for each spin structure, although the partition functions corresponding 
to spin structures of Arf invariant 1 are 0. The sum of these fermionic partition functions, 
multiplied by the factor 1/2^, is the bosonic partition function (151]) . 
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There is a d- dimensional chiral free fermion with invertible super- modular functor when 
d is even. The chiral fermionic partition functions are 

(55) c^^5(9)=n(i+^"^^)'' 

n>0 

(56) c^ii^i(g)=n(i-^""'^)'' 

n>0 

(57) ,Z^'iq) = 2'/'l[il + q-r, 

n>0 

(58) ,Z^{q) = 

(We use the symbol cZ to distinguish the chiral partition functions from the fermionic par- 
tition functions at both chiralities, which are written below in ( 159|) . ( 160|) . (16T]) . (162|) . since 
there is only one label for each spin structure.) The conformal field theory assembled from 
this super-RCFT is isomorphic to the conformal field theory of the bosonised d- dimensional 
fermion described above. 

For d odd, there is no d- dimensional chiral fermion super- conformal field theory as defined 
above. It was proved in [22] that there is no invertible super-modular functor of central charge 
1/2. However, P.Deligne found an even more involved axiomatization using the super Brauer 
group and central simple algebras which does allow a 1-dimensional chiral fermion model. 
Roughly, the super-Brauer group of C is Z/2, and then non-zero element is represented by 
the super-division algebra C = C[a]/ (a^ — 1) where a is odd. Then the value of the modular 
functor (and the Hilbert spaces) are taken in the category of modules over super-division 
algebras, and identification is allowed under super-Morita equivalence. When d is odd, 1- 
dimensional C-modules occur as values of the modular functor when Ramond boundary 
components arise. 

When we put both chiralities together, there is, again, a uniform formula for the fermionic 
(spin-structure dependent) partition functions of the dimensional fermion, regardless of 
the parity of d: 

(59) Z^siQ) = 11(1 + q-^^ni + r^-^r, 

(60) Z^liq) = llil-q-^'^ni-ir-'h', 

n>0 



Z^'{q) = 2''ll{l + q-r{l+rr, 

n>0 

(62) Z^iq) = 

In fact, as a fermionic conformal field theory (i.e. on spin structure at a time), the d- 
dimensional fermion is a tensor product of d copies of the 1-dimensional fermion (and its 
rigorous mathematical construction in fact proceeds in this fashion). 
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(61) 



For d odd, the Hermitian form 

C^C 

has 

(1, 1) = (a, a) = 1, 

which causes, on the R label part of the modular functor, which is given by Deligne's prescrip- 
tion in the form of a (1-dimensional) graded C-module, the additional factor 2, explaining 
the odd power of 2 present in ( 16T|1 . 

Comments and questions about graphs: Physics [8j predicts that there should exist 
a limit formula from the critical dimer model to the 2-dimensional free fermion conformal 
field theory, i.e. the theory which has a bosonic description with modular functor with set of 
labels A2, as well as a fully fermionic description with an invertible super-modular functor. 

On the other hand, physics also predicts [29] that there should exist a limit formula from 
the critical Ising model to the bosonic 1-dimensional fermion theory with the set of labels 
A. We saw that both conformal field theories are closely related yet substantially different. 

The factorization formulas f l45p and f lM|) . and their higher genus analogues, suggest com- 
pelling analogies with decompositions of partition functions of the discrete models into sum- 
mands. This is particularly striking in view of the expression of the dimer and Ising model 
partition functions as a sum of graphs drawn in a suitable way on genus g Riemann surfaces 
as sums of 4^ Pfaffians. In fact, an important point of this paper is that in Sections [2] and 
|3l we have modeled discrete gluing analogues for the fermionic form of the conformal field 
theory, i.e. ([59]), (EOD, ([6l]), (ES]). 

It is not clear whether (or in what situations) the estimates for number of summands given 
by the limit theory are precise: In the conjectured limit, in the fermionic description, fl62p . 
and more generally, the partition functions corresponding to spin structures of Arf invariant 
1 , always vanish (see also [S] ) . This may appear to contradict the result of Norin [27] showing 
that for g = 1, a 3-Pfaffian graph is always 1-Pfaffian, while the limit formula predicts 3 
non-vanishing factors. However, Norin's result |27] is stated in a form not involving weights 
of edges. When translated to our formulation with weights, cancellations may occur, and 
some terms may vanish. This occurs specifically in the case of critical embeddings, even 
without taking the limit. 

Note that in the bosonic description, the interpretation of the limit is somewhat different. 
In this case, the summands are not Pfaffians, but suitable averages of Pfaffians using the 
GSO projection. In this case, in the even-dimensional case, none of the summands actually 
vanish, but subsets of different summands have identical partition functions: in the discrete 
case, then, the question of vanishing is replaced by other linear dependencies. 

In the odd-dimensional case (which is relevant to the Ising model), the situation is more 
interesting: In the bosonic description, one doesn't combine sectors with the same partition 
function, but instead the sectors themselves combine to one sector. Is there a definition of 
corresponding "GSO-projected" summands purely in terms of graphs? 

Alternately, in the fermionic description, is there a graph-theoretic characterization of 
Pfaffians corresponding to the modules over the Deligne super-algebra C? While an ex- 
pression of the Ising partition function as a sum of 4^ Pfaffians is obtained in [21], this 
description uses a modifying construction on the underlying graph. Is there a formula in 
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terms purely of Pfaffian-like summands on the original graph, suitably generalized using 
C-modules discretely mimicking Deligne's theory? 
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